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MIHIIIPIHTIAL 


FOREWORD 


This  report  was  prepared  by  the  Wind  Tunnel  Branch,  Aircraft 
Laboratory,  Directorate  of  Laboratories,  Wright  Air  Development 
Center,  Wright-Patterson  Air  Force  Base,  Ohio.  The  project  was 
administered  under  Research  and  Development  Order  Number  458-429, 
(UNCLASSIFIED),  “Two-Dimensional  and  Axially  Symmetric  Transonic 
Flow,*  with  Gottfried  Guderley  acting  as  project  engineer. 

This  entire  report  is  classified  CONFIDENTIAL  because  it 
contains  informative  and  descriptive  materials  relating  to  the  design 
of  a new  type  of  test-section  wall  for  circular  wind  tunnels.  The 
security  classification  for  such  material  has  been  established  by  the 
National  Advisory  Committee  for  Aeronautics. 
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ABSTRACT 


On  the  basis  of  a simplified  boundary  condition  which  has  been 
derived  in  a previous  report,  the  wind-tunnel  corrections  for  walls 
with  slots  of  different  widths  have  been  determined  by  a method 
analogous  to  that  of  I.  Lotz.  Formulae  have  been  given  for  the  source, 
the  doublet  and  the  infinitesimally  small  horseshoe  vortex.  The  correc- 
tion velocities  at  the  wind-tunnel  axis,  the  average  velocities  at  the 
wind-tunnel  walls  and  the  height  up  to  which  the  jet  will  penetrate  into 
the  slots  have  been  computed  and  are  represented  in  curves. 


PUBLICATION  REVIEW 
This  report  has  been  reviewed  and  is  approved. 
FOR  THE  COMMANDER: 
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Colonel,  USAF 
Chief,  Aircraft  Laboratory 
Directorate  of  Laboratories 


WADC  TR  54-22 


iii 


f 


- COrffw^n  ■ 

table  of  contents 


Section 

I 

n 

HI 

IV 


Expressions  for  the  Singularities  in  Free  Air 

The  Wall  Influence * 

Specialized  Formulae 

Survey  of  th®  Figures  

List  of  References 

Appendix 


Page 

1 

6 

9 

15 

17 

32 


WADC  TR  54-22 


iv 


LIST  OF  ILLUSTRATIONS 


Figure  Page 


1 Velocity  Changes  at  the  Wind  Tunnel  Axis  for  a 

Source  (Functions  Fg  (x)  Versus  x/rD  for  Various 
Values  of  c or  of  C*) . 18 

2 Velocity  Changes  at  the  Wind  Tunnel  Axis  for,  a 
Doublet  (Function  FD  Versus  x/ r for  Various  Values 

of  c or  of  C*)  19 

3 Down  wash  Due  to  the  Wall  Influence  for  a Horseshoe 
Vortex  (Function  Fl  Versus  x/r0  for  Various  Values 

of  c or  of  C*) 20 

4 Average  Wall  Velocity  for  a Source  (Function  Gg 

Versus  Vro  *or  Various  Values  of  c or  of  C*)  21 

5 Average  Wall  Velocity  for  a Doublet  (Function  Gj) 

Versus  x/r0  for  Various  Values  of  c or  of  C*) 22 

6 Average  Wall  Velocity  for  a Horseshoe  Vortex 
(Function  GL  Versus  x/ro  l°r  Various  Values  of  c or 

of  C*) 23 


7 Penetration  of  the  Jet  into  the  Slots  for  Source 

(Function  Kg  Versus  x/ r0  for  One  Value  of  c and  the 
Corresponding  Value  of  C*).  (The  value  C*  is  chosen 
in  such  a manner  that  the  correction  velocity  for  the 


doublet  is  close  to  zero.) 24 

8 Penetration  of  the  Jet  into  the  Slots  for  Doublet 
(Function  Kjv  Versus  x/rQ  for  the  Values  of  c or  C* 

of  Figure  7). 25 

9 Penetration  of  the  Jet  into  the  Slots  for  Lift 
(Function  Kl  Versus  x/rQ  for  the  Values  of  c or  C* 

of  Figure  7) 26 


WADC  TR  54-22 


v 


wiviDlifnM 

LIST  OF  ILLUSTRATIONS  (Continued) 


Figure  Page 

10  Correction  Velocities  at  the  Wind  Tunnel  Axis  for 
x = 0 Versus  C*"*  (Function  Fjq  (0)  Versus  C*“*). 

(The  scattered  points  correspond  to  the  results  of 

the  computation  in  Reference  3 .) 27 


11 

12 


13 

14 


Correction  Velocities  as  in  Figure  10  Versus 

Opening  Ratio  c for  Wind  Tunnel  with  Eight  Slots 28 

Ratio  of  Average  Wall  Velocity  Deviation  to  the 
Velocity  Perturbations  Due  to  the  Wall  at  the  Wind 
Tunnel  Axis  Versus  C* 29 


Ratio  of  Maximum  to  Average  Wall  Velocity 
Deviation  (R  Versus  c) 


Curve  for  the  Determination  of  the  Camber 


Correction 


id  FL 


fo-ry.  =0)  Versus  C* 


30 


31 


NOTE:  Wherever  the  number  of  slots  in  the  wind  tunnel  is  used 
(Figures  1 - 14,  inclusive),  it  is  assumed  to  be  eight. 


WADC  TR  54-22 


vi 


LIST  OF  SYMBOLS 


x,  r,  (o 
x,  y,  z 


n 

c 

C* 

ro 

rl 

F 

CL 

U 

V 

9 

r 

s 


?. 


% 


system  of  cylindrical  coordinates 

system  of  Cartesian  coordinates,  the  z axis  lies  in  the 
plane  cO  = 0 

number  of  slots  in  the  tunnel 

ratio  of  open  to  the  entire  jet  surface  (opening  ratio) 

* - £ ^ ¥ 

radius  of  the  tunnel 

radius  of  half  body  due  to  a source 

wing  area 

lift  coefficient 

free  stream  velocity 

Volume  of  slender  body 

density  erf  the  undisturbed  flow 

circulation 

span  of  horseshoe  vortex 

potential  which  gives  the  deviation  of  the  flow  over  the 
body  in  free  air  from  a parallel  flow 

potential  which  describes  the  wall  influence 

the  potential  which  would  describe  the  flow  over  the 
body  in  the  tunnel,  if  the  simplified  boundary  condition 
(1)  were  valid 


WADC  TR  54-22 


vii 


LIST  OF  SYMBOLS  (Continued) 

^ variable  of  integration  which  occurs  in  a Fourier 

integral 

ag,  aj y aL  functions  which  occur  in  Fourier  integral!  and  which 

are  related  to  the  source,  the  doublet  or  the  horseshoe 
vortex 

FS , Fd,  Fl  functions  which  express  the  wall  influence  along  the 
wind-tunnel  axis  for  the  source,  the  doublet  and  the 
horseshoe  vortex 

Gg,  Gp,  G^  functions  which  express  the  average  wall  velocities  for 
the  source,  the  doublet  and  the  horseshoe  vortex 

HS>  Hd,  Hl  the  height  up  to  which  the  jet  penetrates  into  the  slots 
for  the  source,  the  doublet  and  the  horseshoe  vortex 

KS,  Kjj,  Kl  functions  by  which  the  H's  are  expressed 

R the  ratio  of  the  maximum  of  the  additional  wall  velocity 

to  the  average  wall  velocity  at  x = 0. 


INTRODUCTION 


The  present  report  gives  subsonic  wall  corrections  for  circular 
wind  tunnels  with  longitudinal  slots.  The  most  important  point  is  the 
introduction  of  the  boundary  conditions.  The  simplifications  which  one 
can  make  are  studied  in  a separate  report  (Reference  1)„  There  the 
conclusion  is  reached  that  for  slots  with  a constant  width  which  extend 
from  negative  to  positive  infinity,  the  wall  acts  in  the  average  as  if 
everywhere  the  condition 


with 


<fi-C  rc  f>r  = 0 
C = -Jr  £co  Sun, 


(1) 


is  imposed.  Here  <j>  is  the  potential  which  would  describe  the  devi- 
ation of  the  flow  from  a parallel  flow,  if  Equation  (1)  were  the  actual 
boundary  condition  at  the  wall  (not  only  an  average  expression),  n is 
the  number  of  slots,  r_  is  the  wind-tunnel  radius,  c is  the  ratio  of  open 
area  to  the  entire  surface  of  the  jet  (opening  ratio).  The  above  boundary 
condition  as  an  average  condition  does  not  show  a difference  between 
the  solid  part  of  the  wall  and  the  slot.  So  one  cannot  expect  that  one 
obtains  f>r  =0  along  the  solid  part  of  the  wall  and  <f  =0  in  the 
slots  as  the  actual  boundary  conditions  would  require.  Incidentally  the 
quantity  C*  used  here  is  denoted  in  Reference  1 by  K. 

With  this  boundary  condition  the  wind-tunnel  corrections  can  be 
determined  by  a procedure  which  is  quite  analogous  to  that  of  I.  Lotz 
(Reference  2). 

As  usual  the  body  is  replaced  by  singularities  at  the  tunnel  axis. 
The  following  singularities  will  be  treated: 

1.  the  doublet 

2.  the  source 

3.  the  infinitesimally  small  horseshoe  vortex. 
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The  opening  ratio  c will  be  varied  from  zero  (closed  tunnel)  to  one 
(open  tunnel).  For  these  singularities  the  following  data  will  be  given: 

1.  the  additional  velocities  at  the  tunnel  axis  due  to  the  wall 
influence,  or  the  downwash  at  the  axis 

2.  the  average  values  of  % at  the  wind-tunnel  wall 

3.  the  height  up  to  which,  under  the  present  assumptions,  the 
slots  will  be  filled  with  air  from  the  tunnel. 

The  data  will  be  described  in  a separate  section. 

For  the  computation  the  potential  p is  split  up  into  two  parts -- 
the  first  one,  <pt  , describes  the  flow  in  free  air;  the  second  one,  P 2> 
gives  the  influence  of  the  walls. 

Attention  is  drawn  to  a report  by  Don  D.  Davis,  Jr.  and  Dewey 
Moore,  Reference  6,  which  appeared  after  the  present  report  was 
finished  and  treats  the  same  subject  in  a very  similar  manner. 
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SECTION  I 


EXPRESSIONS  FOR  THE  SINGULARITIES  IN  FREE  AIR 


In  the  following  computations  only  an  incompressible  flow  need 
be  considered  since  the  results  can  readily  be  transformed  to  subsonic 
flow  fields  by  means  of  a Prandtl-Glauert  correction. 

In  this  paragraph  representations  of  the  various  singularities 
which  are  assumed  at  the  axis  will  be  derived  in  the  form  of  Fourier 
integrals.  Once  these  expressions  are  known  the  correction  potential 
can  be  found  nearly  immediately. 

Let  x,  r,  u)  be  a system  of  cylindrical  coordinates  whose 
axis  coincides  with  the  tunnel  axis 

U the  'velocity  of  the  undisturbed  flow 

p the  density  of  the  undisturbed  flow. 

Occasionally  a system  of  Cartesian  coordinates  x,  y,  z will  be  used. 
The  z axis  will  lie  in  the  plane  to  = 0. 

For  an  axisymmetric  flow  the  potential  equation  ^ = 0 has 

the  form 


A,  + 7 ' yxx 


o 


(2) 


We  shall  use  particular  solution  of  this  differential  equation  in  the 
form 


Cp  - Z A r ) c<h  A X 


(3) 
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Here  Z0  denotes  a linear  combination  of  Bessel  functions  of  zero 
order,  A is  a real  quantity.  * Among  these  solutions  the  following 
ones  are  of  particular  importance 


c n"‘  (oXr)  cn*X 


(4a) 


and 


J / i,X  rj  Cc rs 

°'  (4b) 

The  first  functions  vanish  exponentially  as  r goes  to  infinity,  functions 
of  this  type  will  be  used  for  the  representation  of  the  singularities  at 
the  axis.  The  other  particular  solution  is  the  only  one  which  is  regular 
for  r = 0,  it  will  be  used  for  the  representation  of  the  wall  influence. 

To  discuss  the  expression  (4a),  we  determine  the  behavior  of  the 
Hankel  functions  for  small  values  of  r.  With  the  known  formulae  for 
Bessel  functions  (cf.,  e.g.,  Reference  4)  one  obtains 

i H‘"(ar)  ~ -M  (Hr)  -v  - ^ A.  r 

Close  to  the  axis  one  hence  obtains 

' d fi  - - i-7 

dr 


Thus  it  follows  that  the  expression  (4a)  gives  a flow  field  with  a 
source  distribution  along  the  x axis,  for  which  the  source  strength 
per  unit  length  is  given  by 

- ^ q c*s  -A  Jf 
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A suitable  superposition  of  the  expression  (4a)  will  give  a source 
of  unit  strength  at  the  origin.  The  Fourier  integral  for  a <T  -function 
of  weight  one  at  the  origin  is  given  by 


Thus  one  obtains  as  the  desired  expression  for  the  potential  of  a source 
with  unit  strength  at  the  origin 


X 


* (5) 


In  the  notation  of  <f>  the  subscript  S indicates  that  the  expression 
refers  to  a source;  later  the  subscripts  D for  the  doublet  and  L for  the 
lifting  element  will  be  introduced.  The  subscript  one  with  Q>  indicates 
that  the  expression  refers  to  free-stream  conditions.  The  last  expres- 
sion converges  everywhere  except  at  the  axis;  actually,  it  will  be  used 
only  at  the  wind-tunnel  wall.  The  corresponding  closed  form  for  the 
potential  would  be 


/ t 

vW+W ' 


The  potential  for  a doublet  is  obtained  by  a differentiation  with 
respect  to  x and  a change  of  sign.  Then  the  moment  of  this  doublet  has 
the  value  one.  Thus  one  obtains 


OO 


(6) 


According  to  a well-known  result  of  slender-body  theory  the 
doublet  strength  equivalent  to  a body  of  volume  V is  given  by  -XIV  p 
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For  the  sphere  one  obtains  a very  similar  result,  namely,  -3/2*UV^  . 
For  intermediate  bodies  one  will  use  a correction  factor  with  the 
expression  -UV  q which  depends  upon  the  body  shape  (see  e.g., 

Reference  5,  Figure  2). 

To  arrive  at  the  potential  of  an  infinitesimally  small  horseshoe 
vortex  with  a distance  of  the  two  trailing  vortices  s and  a circulation 
r , we  investigate  the  behavior  of  its  potential  close  to  the  axis. 

The  horseshoe  vortex  may  lie  in  the  plane  z = 0.  If  s is  small,  then  at 
a station  x sufficiently  away  from  the  airfoil  one  will  have  essentially 
the  two  dimensional  potential,  i.e., 

£ Jose -arc ) 


For  a small  value  of  s this  gives 


£ £ 
2X  r* 


The  result  indicates  that  along  the  positive  x axis  the  potential  of  a 
horseshoe  vortex  behaves  like 

S f*  .JL  cot  to 

Zx  r 


Along  the  negative  x axis  it  is  obviously  regular.  To  obtain  an 
expression  for  such  a potential  we  first  decompose  it  into  a symmetri- 
cal and  into  an  antisymmetrical  part,  with  respect  to  the  plane  x = 0, 
i.e.,  into  the  potential  of  two  vortices  of  opposite  sign  with  strength 

£ which  extend  from  negative  infinity  and  into  the  potential  of  a 

vortex  system  which  agrees  with  the  symmetrical  part  for  positive 
values  of  x and  has  the  opposite  sign  for  negative  values  of  x.  The 
potential  of  the  first  part  is  given  by 

. -L  cos  <o 
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The  potential  of  the  second  part  will  be  obtained  by  a super -position 
of  expressions  of  the  form 

- H * \ y ) COS  Co  • SlA*s  fx 


They  are  solutions  of  the  potential  equation.  Along  the  x axis  one  has 


5 


CO 


The  Fourier  integral  for  the  unit  step  function  is 


OO 


Thus  one  obtains  for  the  potential  of  the  antisymmetric  vortex 
distribution 


OO 


•srctuo  f(_  y°>{a  rj.  n ij  yj 

4 X r0  J ' yc 


The  lift  produced  by  the  bound  vortex  is 

S Tp  U 
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Thus  the  potential  for  an  infinitesimally  small  horseshoe  vortex  with 
unit  lift  is  given  by 


Cos  <o 

f lt4zr0 


( 


5 


O 


An  alternative  procedure  to  derive  this  expression  is  the  follow- 
ing: one  starts  with  an  even  source  distribution  in  that  region  of  the 
plane  z = 0 which  is  enclosed  by  the  horseshoe  vortex.  The  source 
strength  per  unit  of  area  may  be  one.  The  potential  of  this  flow  is 
continuous  everywhere,  its  derivative  in  the  z direction  makes  a unit 
jump  along  this  sheet.  Since  the  derivative  of  the  potential  with  respect 
to  z fulfills  also  the  potential  equation,  this  derivative  represents  a new 
potential  function  which  is  continuous  everywhere  except  in  the  portion 
of  the  plane  z = 0 which  is  enclosed  by  the  horseshoe  vortex.  This 
expression  then  represents  the  desired  potential  of  the  horseshoe  vortex. 
Naturally,  this  procedure  can  be  used  for  an  infinitesimally  narrow 
horseshoe  vortex.  The  potential  of  the  source  distribution  is  at  some 
distance  from  the  vortex,  that  of  a line  source  whose  discharge  of  mass 
per  unit  length  is  given  by  s (s  was  introduced  as  the  span  of  the  horse- 
shoe vortex).  This  expression  is  given  by 


The  derivative  with  respect  to  z gives  then  the  potential  of  the  horseshoe 
vortex.  Thus  one  obtains  again  the  previous  formula. 


SECTION  n 

THE  WALL  INFLUENCE 


One  must  now  determine  a solution  <j>j_  of  the  potential  equation 
which  is  regular  inside  of  the  tunnel  and  which,  if  added  to  the 
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expression  for  the  potential  of  a source  a doublet  and  a horseshoe 
vortex,  fulfills  the  boundary  condition  (1)„ 

In  the  following  equations  ag(  ^ )>  ao(  ^ ) and  a^(  ^ ) will  be 
suitable  functions  of  A , used  to  express  the  potential  by  means  of  a 
Fourier  integral. 

The  correction  potential  for  a source  may  be  expressed  by 


tu  * A <A)  Jcfc}  rj  **  to  Vdx 

¥ * (ID 

The  function  ag(  \ ) is  determined  by  inserting  the  expression 

into  the  boundary  condition  (1).  <p/s  is  given  by  Equation 
(5).  At  the  txmnel  boundary  one  has  r = rQ.  One  thus  obtains 


n ■ c'l (c  fax 

Vr°  X r~ra 


- % 


(A)  (4  (a)  / IS  C*d  (Jo  fix  $jJ  cn  (A%)JA  .0 


Hence  one  finds 


i H 7‘V  + Cr,  xiritr0  (a 

9 - ^ 


CL-Cdi)  , 

/ C-r.i  Re(c'H> 


'dr 
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Using  the  formula  for  the  derivative  of  a Bessel  function  and  inserting 
the  resulting  expression  for  ag(  A ) into  Equation  (11)  one  obtains 


hr 


4fr„  X 


The  corresponding  egression  for  the  doublet  can  be  found  in  an 
analogous  manner,  but  it  can  be  found  more  quickly,  if  we  remember 
that  the  potential  is  the  negative  derivative  of  s with  respect 
to  x.  The  boundary  condition  is  prescribed  along  a surface  r = r0  = 
constant;  therefore,  it  can  be  differentiated  with  respect  to  x.  Accord- 
ingly, if  one  differentiates  the  entire  expression  for  the  source 

+ s with  resPect  to  x> the  result  will  fulfill  the  boundary 
condition  (1).  Thus 

A* . * ~ ^ 


or  explicitly 


H roX;C  J 


lHqUI)  - c*J( 


. J (tX  Sth (x\  )dx 


(13) 


For  the  lift  case  the  investigation  is  quite  similar.  We  set 


. cen 

fit-  ~ fUtZ’o 


tit f feuf-WH.) ) 
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where  / is  a constant  which  must  be  determined.  Inserting  the 
expression  <b  + CP. , into  the  boundary  conditions  one  obtains 
finally  'L‘  ' 

9 C«a>  fC*-!  r J(,-c')f-Kh)-C\cUV) 

TlL  ' £ ■ !(i-c')(-J,u\)+c^  Jji\) 


SECTION  in 

SPECIALIZED  FORMULAE 

From  the  expressions  (5),  (6),  (9),  (12),  (13)  and  (14)  one  obtains 
without  difficulty  the  quantities  which  are  of  practical  interest. 

Since  J (o)  = / 

one  finds  for  the  wind-tunnel  axis  r = 0 


C)J(  ' 


/ 


o JQ  fcX)  -f.  c* X f-iJ,  fa) 


• W /Cl  jz)  cf  X 
(15) 


(16) 
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For  reasons  of  symmetry  the  horseshoe  vortex  does  not  produce 
additional  values  of  fc*  along  the  wind-tunnel  axis,  but  it  will  give  a 
downwash.  The  downwash  is  determined  as 

d2* 

cf  HL 


Now  one  has 


. i 

r 'doo 


• S(/hA) 


cfz  *dr 

The  series  development  of  Bessel  functions  yields  for  small  values  of  A 
Thus  one  obtains  from  Equation  (14) 


'c*-/ 


j (tc!)L!L (QL  — I a'  /i  J )dX 

1/T'  p U^zt‘  ( c*-n  * J (i-cVHJuv-fC**  TjcX)  r° 

6 (17) 

Equations  (15),  (16)  and  (17)  may  be  written  in  a form  which  puts 
the  geometrical  quantities  into  better  evidence.  As  is  known  a source 
in  a parallel  flow  corresponds  to  the  flow  over  a half  body.  If  the 
radius  of  the  half  body  is  then  the  source  strength  is  Tc  r/p  U . 

Thus  the  expression  (15)  can  be  rewritten  as 


U -£■/$&/ 
0 


along  the  axis 


(18) 


with 


¥ ]j  „■»,/  T,.< j • -z)  d A 


J„ux  +c*x(-iJ/a>) 


(18a) 
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c^l-t 


It  was  mentioned  previously  that  a slender  body  of  volume  V 
is  equivalent  to  a doublet  of  strength  -UV  ^ . Hence  one  obtains  from 
Equation  (16) 


along  the  axis 


(19) 


with 


^ 7 i H0°ax)  - c*)/- Hf'h'to) 
' ' *■ y j.  ea ) + c 


■tfctsfi  f Jett 

ro 

(19a) 


Finally  one  can  express  the  lift  L by  means  of  the  lift  coefficient  c^ 
and  the  area  of  the  model  F 

L = F*cl.1/2  (o  U2 

Then  one  obtains  from  Equation  (17)  for  the  change  of  the  angle  of 
attack 


« 


Fcl 

<nr> 


along  the  axis 


(20) 


with 


t [*J ' z Tfw  7 F5  ^7 ’ A 

0 (20a) 


The  average  velocities  at  the  walls  are  obtained  by  putting  r = rQ  in 
the  expression 


h - /?*  * a* 


Thus  one  finds 
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& - * f.W«  - 


By  the  formulae  on  page  144  of  Reference  4 the  last  expression  can  be 
simplified  and  one  obtains 


%s  **) 


along  the  wall 


(21) 


where 


Co 


aexj.  £ / 


Jt  • •«>*$  f)  J'X 

'o 


(21a) 


similarly 


% 


where 


“vf  s^> 


along  the  wall 


(22) 


OO 


£ 

J) 


C 


***■;  J0  + ri  j-cJavJ 


J • cos,  (J  ) d ^ 


(22a) 


and 
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along  the  wall 


(p  = ^ 7— — (T  fa  / 

Axz  /r  r/  z- 

where 

^ _ 2jT  * f ^ ^ j * 

L ' * 1 (l-c'K-J'/m+C^X- • 

(23a) 


Furthermore,  the  height  up  to  which  the  wind-tunnel  air  enters 
the  slots  will  be  determined.  It  is  found  by  an  integration  of  the  stream- 
line slopes  in  the  slots.  The  procedure  is  the  following  equations.  From 
Equations  (21)  and  (21a)  one  obtains  by  an  integration 


jf 


- OO 


This  is  the  average  potential  along  the  wall.  Hence  one  obtains  L> 
means  of  the  boundary  condition  (1) 


This  is  the  average  value  of  <f>Y  . For  reasons  of  continuity  the  value 
of  (j\  in  the  slot  is  this  average  value  multiplied  with  the  reciprocal 
of  the  opening  ratio  c.  Thus 


oc 


- Oo 


Let  us  denote  the  heights  up  to  which  the  slots  are  filled  by  Ho,  Hr* 
and  Hl»  One  then  finds  0 
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(24) 


n c £ vJ 


where 


^ U)  * A*  /*** 

L ^ Oo  -°o 


Here  j*  and  are  variables  of  integration. 
Similarly 

. ■ i ({*)•%*> 


where 


^ Cot) 


C 


VP*f 


■ Oo  -Oo 


and 


M,  _ ^ P c. 


rz  c <fzr6*  A 


K (x)  • ^Af<b 


where 


/£  £*;  55 


* A*  /k  OV  A 


- oo  — Oo 


(24a) 


(25) 


(25a) 


(26) 


(26a) 
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SECTION  IV 


SURVEY  OF  THE  FIGURES 


The  functions  necessary  to  determine  the  wall  influence  along 
the  wind-tunnel  axis  are  shown  in  Figures  1,  2 and  3 for  different 
values  of  C*  as  a function  of  x/ro*  Besides  the  values  of  C*  the  values 
of  c for  a wind  tunnel  with  eight  slots  have  been  given.  According  to 
Equations  (18)  and  (19)  Figures  1 and  2 give  the  velocities  along  the 
wind-tunnel  axis  caused  by  the  wall  influence  for  a source  and  a doublet 

Figure  1 shows  that  a slotted  tunnel  does  not  require  a Mach 
number  correction  due  to  the  wake  of  the  body.  Such  a Mach  number 
correction  will  occur  if  the  Mach  number  change  due  to  the  wall  far 
upstream  of  the  model  is  different  from  the  Mach  number  change  at  the 
location  of  the  model.  However,  this  happens  only  for  a closed  tunnel 
and  not  for  slotted  tunnels,  since  the  presence  of  the  slots  enforces  far 
upstream  and  far  downstream  a Mach  number  which  is  determined  by 
the  pressure  outside  of  the  slots. 

In  a closed  tunnel  a correction  of  the  drag  due  to  the  pressure 
gradient  caused  by  the  wake  may  have  some  importance.  Figure  1 
shows  that  in  slotted  tunnels  this  correction  would  be  only  a fraction 
of  the  correction  in  a closed  tunnel. 

Figure  2 gives  the  influence  of  the  thickness  of  a body.  One  is 
primarily  interested  in  the  Mach  number  correction  at  the  location  of 
the  model.  Therefore,  Figure  10  shows  the  function  Fj)  versus  1/C*. 
The  scattered  points  are  the  points  which  can  be  found  in  the  corres- 
ponding figure  of  Reference  3.  In  Figure  11  the  same  quantity  has  been 
plotted  over  the  opening  ratio  c for  a wind  tunnel  with  eight  slots.  This 
graph  shows  very  distinctly  that  even  an  extremely  small  opening 
ratio  gives  a very  appreciable  deviation  from  the  values  of  a closed 
tunnel.  The  behavior  of  the  curve  11  at  c = 0 is  studied  in  an  appendix. 
Practically,  non-linear  effects  would  give  a more  gradual  transition 
for  small  values  of  the  opening  ratio. 

Figure  3 gives  essentially  the  downwash  along  the  axis  except  for 
the  quantity  (C*  - 1)/(C*  + 1)  in  Equation  (20).  At  the  location  of  the 
model  the  downwash  is  determined  by  this  quantity  only.  Therefore, 
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Figure  3 illustrates  only  how  the  downwash  changes.  .Eventually  the 
results  could  be  expressed  as  a camber  correction.  For  this  purpose 
Figure  14  shows  for  x = 0 


dlL 


versus  1/C*. 


It  is  remarkable  that  the  value  of  C*  for  zero  Mach  number  correction 
(Figure  10)  and  for  zero  camber  correction  (Figure  14)  are  not  too 
different;  thus  camber  correction  and  Mach  number  correction  can 
practically  be  eliminated  simultaneously. 

The  average  wall  velocities  can  be  obtained  from  the  graphs  4,  5 
and  6 by  Equations  (21),  (22)  and  (23).  The  order  of  magnitude  of  these 
perturbations  is  quite  appreciable  and  suggests  caution  in  the  applica- 
tion of  the  present  results  in  the  transonic  range.  To  show  the  relative 
orders  of  magnitudes,  the  ratio  of  the  Mach  number  correction  along 
the  tunnel  axis  at  the  location  of  the  model  to  the  Mach  number  devia- 
tion at  the  wall  from  the  free-stream  Mach  number  (also  for  x = 0)  has 
been  plotted  for  a doublet  in  Figure  12  versus  C*.  In  the  small  range 
of  values  C*  shown  there,  this  ratio  is  quite  small. 


The  average  wall  velocity  which  is  obtained  from  Equations 
(21),  (22)  and  (23)  cannot  be  measured  directly.  However,  one  can 
measure  the  maximum  wall  velocity  by  pressure  holes  in  the  middle  of 
the  solid  part  of  the  wall.  The  ratio  between  the  average  and  the  maxi- 
mum velocity  deviation,  may  be  denoted  by  R,  depends  upon  the  opening 
ratio  c and  can  be  found  from  the  consideration  of  Reference  1.  This 
ratio  is  represented  in  Figure  13. 

The  penetration  of  the  wind-tunnel  air  into  the  slots  can  be  found 
from  the  curves  7,  8 and  9 by  means  of  the  formulae  (24),  (25)  and 
(26).  Reference  1 showed  that  the  penetration  of  the  tunnel  air  into  the 
slot  may  have  a very  essential  influence  with  respect  to  the  boundary 
condition.  These  graphs  can  be  used  in  order  to  obtain  an  insight  into 
possible  errors  for  a given  case. 

In  general,  it  is  remarkable  that  a relatively  small  opening  ratio 
suffices  to  make  the  behavior  of  the  slotted  tunnel  quite  similar  to  that 
of  an  open  one.  For  values  of  C*  for  which  the  Mach  number  correc- 
tions at  the  axis  are  small,  the  results  are  not  too  sensitive  against  a 
change  of  the  constant  C*.  This  is  quite  encouraging  from  the  practical 
point  of  view. 
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Figure  2:  Velocity  Changes  at  the  Wind  Tunnel  Axis  for  a Doublet 

(Function  FD  Versus  x/r0  for  Various  Values  of  c or  of  C*) 
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Figure  3:  Downwash  Due  to  the  Wall  Influence  for  a Horseshoe  Vortex 
(Function  Fl  Versus  x/r0  for  Various  Values  of  c or  of  C*) 
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Figure  6:  Average  Wall  Velocity  tor  a Horseshoe  Vortex  (Function 
Versus  rj  rQ  for  Various  Values  of  c or  of  C*) 
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Figure  7:  Penetration  of  the  fet  Into  the  Slots  for  Source  (Function 
Kg  Versus  x/  r^for  One  Value  of  c and  the  Corresponding 
Value  of  C*).  (The  value  C*  is  chosen  in  such  a manner 
that  the  correction  velocity  for  the  doublet  is  close  to  zero.) 
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Figure  8:  Penetration  of  the  3&t  into  the  Slots  for  Doublet  (Function 
Kp  Versus  x/ rQ  for  the  Values  of  c or  C*  of  Giure  7) 
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tration  of  the  Jfet  into  the  Slots  for  Lift  (Function  KL 
us  x/rQ  for  the  Values  of  c or  C*  of  figure  7) 
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APPENDIX 

THE  BEHAVIOR  OF  A CERTAIN  INTEGRAL  FOR  LARGE  VALUES  OF  C* 


To  study  the  integral 


fd\ 


we  first  form  the  difference  between  its  value  for  an  arbitrary  value 
of  C*  and  its  limiting  value  for  C*  * oo.  One  obtains 

fcHj'hxi-c'xf-H'V*))  f f-  £<n 

1 JC(C\)  -tc'Xf-cJM)  / 

, fcdo'nXll-i 

°(U  \)+c*x  HJ,cai)  f-cj,(cxi) 

_ j.  ~ 3 JO, 

* / UfiX) +c*x(-£. 

Now  we  introduce 


Then  the  above  expression  is  given  by 

. fJLAAl — 

y A 2 > r * 1 j j f-  i 7 // 


!)jf- 4 U ) 
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Its  derivative  with  respect  to 


/ 


is 


i 7 J d * i x <n 


To  evaluate  the  expression  for  f - 0 is  not  possible  immediately, 
since  then  the  integrals  will  not  converge.  Therefore,  we  split  the 

-/* 

O 

which  extend  from  £ to  infinity  it  is  allowed  to  put  f = 0.  In  the 
other  integral  the  Bessel  functions  can  be  developed  with  respect  to 
X and  one  obtains 


integrals  up  into  two  parts:  J* 


/ 


In  the  integrals 


i f tjXdl  A / 1/>  x d\ 

*o  QWm  " 


% 

H 0 then  the  upper  limit  of  the  integrals  in  terms  of 

as  variable  of  integration  tends  to  infinity.  The  first  integral  gives 
7£~/2,  the  second  one  & /4.  Thus  the  derivative  of  the  above  expres- 
sion with  respect  to  J>  is  2 • YT. 


Accordingly,  the  curve  of  Figure  10  behaves  at  1/C*  0 as 

constant  + constant  C*-*'  . The  curve  in  Figure  11  behaves  for  c -t  0 
as  constant  + constant  (log  c)"V^, 
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